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Measurements of mass flux are obtained in a vitiated supersonic ground-test facility using a sensor based on line- 
of-sight diode laser absorption of water vapor. Mass flux is determined from the product of measured velocity and 
density. The relative Doppler shift of an absorption transition for beams directed upstream and downstream in the 
flow is used to measure velocity. Temperature is determined from the ratio of absorption signals of two transitions 
(X.j = 1349 nm and X 2 = 1341.5 nm) and is coupled with a facility pressure measurement to obtain density. The 
sensor exploits wavelength-modulation spectroscopy with second-harmonic detection for large signal-to-noise ratios 
and normalization with the If signal for rejection of non-absorption-related transmission fluctuations. The sensor 
line of sight is translated both vertically and horizontally across the test section for spatially resolved measurements. 
Time-resolved measurements of mass flux are used to assess the stability of flow conditions produced by the facility. 
Measurements of mass flux are within 1.5 % of the value obtained using a facility predictive code. The distortion of the 
wavelength-modulation spectroscopy lineshape caused by boundary layers along the laser line of sight is examined 
and the subsequent effect on the measured velocity is discussed. A method for correcting measured velocities for flow 
nonuniformities is introduced and application of this correction brings measured velocities within 4 m/s of the 
predicted value in a 1630 m/s flow. 


Nomenclature 

a = frequency-modulation depth, crrr 1 

c = speed of light in vacuum, m/s 

E" = lower state energy of absorption transition, cm 

f = modulation frequency, Hz 

I a = average laser intensity 

I v = spectral intensity, W/cm 2 s _1 

i a = linear intensity modulation amplitude 

k v = spectral absorption coefficient 

L = path length, cm 

m = modulation index 

P = pressure, atm 

S = linestrength, ever 2 / atm 

T = temperature, K 

U = velocity, m/s 

x = spatial coordinate along line of sight, cm 
= spectral absorbance 
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Av = frequency shift, cm -1 

6 = crossing half-angle between two beams for Doppler shift 

measurement 

v = optical frequency, cm~‘ 

v = average laser frequency, enr 1 

v 0 = transition linecenter frequency, cm~' 
p = density, kg/m 3 

f = lineshape function 

i/f| = linear intensity-frequency phase shift, rad 

I. Introduction 

A CCURATE measurements of gas dynamic properties such as 
mass flux, temperature, and velocity are necessary for the 
characterization of ground-test facilities and for quantifying perfor- 
mance of advanced propulsion systems. In particular, mass flux is an 
important component in calculations of engine parameters such as 
thrust and specific impulse for airbreathing propulsion systems; in a 
supersonic combustion ramjet (scramjet) for example, a measure- 
ment of mass capture to within 1% uncertainty is necessary for a 1% 
uncertainty in the specific impulse [1], Hence the development of 
accurate and robust sensors for the measurement of mass flux is 
critical for both in-flight and ground-testing of propulsion systems; 
this work focuses on the development of such a sensor for vitiated 
ground-test facilities, which are often used to reproduce hypersonic 
flight conditions, either by matching flight enthalpy or velocity. 
A common method for generating high-enthalpy flow is by 
combustion-heating the air, then expanding the flow through a 
supersonic nozzle with optional 0 2 replenishment [2], Conventional 
invasive mass-flux probes tend to disturb the flow and generate shock 
structures and may have limited ability to survive long-duration runs 
[3] . Noninvasive optical sensors are an ideal solution for testing in this 
type of environment; additionally, the high water mole fraction due to 
combustion allows for the use of widely available tunable diode lasers 
to access H 2 0 absorption transitions in the 1.4 pm range. 

Tunable diode laser absorption spectroscopy (TDLAS) has been 
proven as a robust, noninvasive measurement technique for the harsh 
environments commonly experienced in high-speed or combusting 
propulsion flows. Mass-flux sensing via TDLAS was pioneered by 
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Philippe and Hanson [4] and patented for thrust measurements in 
1993 [5], Subsequent TDL mass-flux sensors have been deployed for 
a variety of field tests, including a commercial turbofan (PW6000) 
inlet at Pratt and Whitney [6], a full-scale Pratt & Whitney F-100 
engine in an open ground-test stand at the NASA Dryden Flight 
Research Facility [7], and a model scramjet combustor at Wright- 
Patterson Air Force Base [3]. A TDLAS sensor employing 
wavelength-modulation spectroscopy (WMS) for measurement of 
mass flux using water vapor absorption was previously tested in a 
subsonic tunnel with uniform flow [8]. Optimization of this sensor 
for high-enthalpy flows is discussed here. The sensor was subse- 
quently deployed in a combustion-heated supersonic ground-test 
facility [NASA Langley Direct-Connect Supersonic Combustion 
Test Facility (DCSCTF)] for spatially and temporally resolved mass- 
flux measurements in the tunnel isolator section. These results are 
compared with the facility predictive code and computational fluid 
dynamics (CFD) solutions to assess sensor accuracy and to examine 
the influence of the nonuniformities predicted by the CFD on the 
sensor-determined values. 

TDLAS measurements are based upon Beer’s law, where 
attenuation of the transmitted light is path-integrated along the line of 
sight. When temperature, pressure, and composition are uniform in 
the absorption path, the detected lineshape is a well-defined function 
of the gas properties. The ratio of two lineshapes provides a means to 
determine temperature, and velocity can be inferred from the 
Doppler shift of a lineshape. However, the lineshape can become 
distorted in the presence of nonuniform pressure, temperature, 
velocity, or gas composition. This distortion is a function of the 
degree of nonuniformity and the pressure and temperature depen- 
dence of the selected transition. Distortion of the lineshape can shift 
the apparent transition linecenter, producing an error in the measured 
Doppler shift. Spatially resolved CFD is used here to evaluate the 
effects of nonuniformity on a line-of-sight (LOS) measurement with 
an eventual goal of developing design rules to minimize the influence 
of nonuniformities. This analysis is used to adjust the TDLAS 
measurements to account for nonuniformity in the DCSCTF flow. 

II. Simulation of Path-Integrated WMS Line 
Shapes from Nonuniform Flowfields 

A. Theory 

TDLAS is most accurate under conditions where the absorbing 
medium is uniform along the beam path. Gradients in velocity, 
temperature, pressure and composition, as well as physical obstruc- 
tions and stagnant regions, can cause distortion to the detected 
lineshape. In this section, the equations describing the detected 
lineshape in a nonuniform absorbing medium will be developed. 

As described in [8], the 1 /-normalized WMS-2/ signal (WMS- 
2//1/) enables rejection of non-absorption-related transmission 
fluctuations and also improves velocity resolution. This sensor was 
demonstrated recently in a high-uniformity, low-speed tunnel to 
accurately measure velocities within 0.5 m/s of the facility set point 
[8], In an earlier study, with an 0 2 sensor, Lyle et al. [6] demonstrated 
that at low velocities, the distortion to the detected lineshape caused 
by boundary layers is nearly negligible. Knowing the fraction of the 
beam path passing through stagnant regions of absorption, a simple 
linear correction can be applied to remove the effects of the low- 
velocity regions [6,7], This technique was demonstrated at subsonic 
velocities by simulations assuming simple step-function changes 
in velocity. However, at high velocities, the larger Doppler shift 
results in more noticeable distortion in the lineshape. Actual flows 
possess continuous gradients in velocity, temperature, pressure, and 
mole fraction that are poorly described by a step function. To 
accurately simulate lineshapes produced in nonuniform flows, we 
begin by considering the equation of radiative transfer in the case of 
absorption only: 


d/„(x) = — k v (x)I v (x)dx (1) 

Here, the spectral intensity l v and the spectral absorption coefficient 
k v depend on position in the absorbing medium, a. Next, we integrate 


with respect to a and use the boundary condition /„( 0) = I v a to 
obtain the transmitted intensity I v : 

h = h,o ex P {^j -K d-r j exp(— a„) (2) 

In the case of a uniform absorbing gas sample, Eq. (2) reduces to 
the familiar form of Beer’s law. Instead, we retain the integral from 0 
to L, the path length, recognizing that the spectral absorption 
coefficient may vary along the laser line of sight. The spectral absor- 
bance, a v , introduced here will be examined in more detail below. 
The transmitted intensity along a nonuniform absorption path 
defined in Eq. (2) will be refined as the equations needed to simulate 
the WMS lineshape are developed. 

A brief introduction to WMS theory is now given; nomenclature 
will follow [8-11]. The injection current of the diode laser is driven 
with a high-frequency sinusoid superimposed on a lower-frequency 
wavelength scan. Current modulation results in simultaneous 
modulation of laser frequency and intensity, and the instantaneous 
values are given by 

v(t) = v + a cos(2?r/t) (3) 

foO) =Io + b cos(2jt ft + ft) (4) 

In Eq. (3), v(t) is the laser frequency (cnr 1 ) averaged over the 
modulation, / is the modulation frequency (Hz), and a is the 
modulation depth (cm -1 ). In Eq. (4), /],(/) is the laser intensity 
averaged over the modulation, i 0 is the linear intensity modulation 
amplitude, and is the linear phase shift between intensity and 
frequency. Equations (2-4) can now be combined as follows: 

I,(t) = [/ 0 + *o cos(2jr/t + ft)] exp[— a(v + a cos(2jr ft))] (5) 

Recognizing that the function is even and can be expanded using 
the Fourier cosine series, we obtain 


OO 


exp[— a(v + a cos(2nft))] = 2, Hk( v > a ) cos(2nkt) 

(6) 

H 0 (v, a) = 

k=0 

1 f * 

— J exp[— a(v + a cos(f))] d£ 

(7) 

1 



H k (v, a) = - 

/ exp[— a(v + a cos(£))J cos (kt-) d£ 

(8) 

nj 

—n 



^ S(x ) • P(x) ■ x(x) • (p(x) dx 
0 

(9) 


Note that in Eqs. (7) and (8), the integration variable is designated 
as £ instead of 6 as in [8-111. This is due to the fact that 6 appears in 
later equations as the crossing half-angle used in crossed-beam 
Doppler shift velocimetry. In Eq. (9) the absorbance has been 
expressed as the integral of the product of linestrength S 
(cm -2 atm -1 ), pressure P (atm), mole fraction of the absorbing 
species, /, and lineshape function f v , with each parameter being 
allowed to vary with a, the position along the line of sight. In the case 
of a uniform line of sight, the absorbance in Eq. (9) reduces to k v L, 
and the previous equations are identical to those in [3,6—8] . 

The transmitted intensity then passes through a lock-in amplifier, 
which consists of a mixer and a low-pass filter, producing an X and Y 
component at the output: 


% = 


Gin 


H 2 + j(Ht+H 3 )co S ft 


(10) 


T 2/ = - 


Glo 
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b 
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(Hi — H 3 ) sin ft 
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Here, G, the detector gain, now appears in these equations. The 2/ 
signal is given by the root sum square of the two components: 

S 2 f = y/x%f + Y 2 2f (12) 

As described in [8,9,11-13], normalization with the If signal is 
used to remove the dependence of the 2/ signal on detector gain and 
average laser power. Similar equations for the 1/ signal are 
produced: 



Rif = yjAf + Y u ( 15 ) 

Equations (7-15) provide the basis for simulating the path- 
integrated WMS-2// I f lineshape along a nonuniform line of sight. 
Note that Eqs. (10-15) are the same as those applied to cases with a 
uniform line of sight; only the formulation of the absorbance in 
Eq. (9) is affected by nonuniformities in the beam path. The overall 
nonuniform lineshape includes the effects of different pressures, 
mole fractions, and temperatures (through the linestrength and 
lineshape functions) experienced at each point along the line of sight. 
The analysis will now proceed to consider the effects of Doppler- 
shifting of the lineshape caused by flow velocity. 

A typical setup for TDLAS velocimetry is shown in Fig. L One 
beam is directed upstream in the flow while the second is directed 
downstream. The detected lineshapes on either beam will experience 
frequency shifts of opposite sign, and the axial velocity, U, can be 
extracted by measuring the relative shift between the two lineshapes: 

Av/v D = 2 sin 9 ■ U/c (16) 

Here, Av is the frequency shift (cm -1 ), v a is the unshifted linecenter 
frequency (cm -1 ), 0 is the crossing half-angle seen in Fig. 1, and c is 
the speed of light (m/s). From Eq. (16 ) it can be seen that lineshapes 
along the beam path may experience slightly different frequency 
shifts caused by nonuniformity in the axial velocity, U, along the line 
of sight. To properly simulate the effects of velocity nonuniformity, 
Eqs. (7) and (8) are modified as follows: 


1 f” 

H 0 (v,a )= — / exp[— a(v + Av(x:)/2 + a cos(f))] d£ (17) 


1 f” 

H k (y , a) = — I exp[— a(v + Av(x)/2 + a cos(f))] cos (££) d£ 

X J-n r 

(18) 

The relative frequency shift Av [as defined in Eq. (16)], has now 
been included to indicate that the absorbance may be frequency- 
shifted due to velocity along the beam path (relative frequency shift is 
divided by two, since these equations consider absorption on a single 
beam). The sign of A v is positive for a beam directed downstream 
and negative for a beam directed upstream; the frequency shift is 
again a function of x, indicating that different velocities at points 
along the beam path will cause different magnitudes of frequency 
shift. Equations (9-18) now completely define the path-integrated 
WMS-2// 1/ lineshape allowing for temperature, pressure, mole 
fraction, and velocity nonuniformities along the laser line of sight. As 
described in the next section, these equations will be used to simulate 
WMS lineshapes in a nonuniform flow. 

B. Lineshape Modeling 

A CFD solution for the flow through the DCSCTF isolator where 
the mass-flux sensor was deployed was used to simulate lineshapes 
for lines of sight through nonuniform flow. The solution provides 
temperature, pressure, density, and uvw velocities in three dimen- 
sions. Hence this data can be compared with spatially resolved 
measurements taken by translating the TDLAS LOS during a run. 
Figure 2a shows the sensor mounted to stages to translate the sensor 
LOS. The setup can be configured to translate either vertically or 
horizontally across the duct. 

In the vertical translation configuration of Fig. 2b, the lasers cross 
in a horizontal x-z plane, which was then translated vertically (y axis) 
across the duct. Similarly for the horizontal translation configuration 
of Fig. 2c, the lasers cross in a vertical x-y plane and were translated 
horizontally (z axis) across the duct. Because the sensor collects 
measurements as it translates across the duct, planes of CFD data 
(outlined in black) encompassing all laser lines of sight (indicated by 
the dashed arrows of Figs. 2b and 2c) can be evaluated. Spatial 
resolution was limited by the finite beam size of the laser; hence the 
CFD data, which was produced on a very fine spatial grid, was 
averaged over the beam diameter (approximately 1 mm). The CFD 
solution for pressure along the laser LOS for the vertical translation 
configuration is displayed in Fig. 3a. The planes of data for 



Fig. 1 Two-laser frequency-multiplexed WMS sensor for mass flux at H 2 0 wavelengths 7 , and X , (~1349 and 1341.5 nm, respectively). Light from the 
two lasers is combined into a single fiber and then split to be directed upstream and downstream in the supersonic flow with a crossing angle 29. Velocity is 
determined from the relative Doppler shifts of the absorption lineshape, and gas temperature is determined from the ratio of the two absorption signals. 
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Fig. 2 NASA DCSCTF isolator: a) section with TDLAS mass-flux sensor configured for vertical scanning, b) CFD geometry (symmetry about vertical 
axis is assumed) with vertical translation configuration shown, and c) CFD geometry with horizontal translation configuration shown. 




•m 1600 

j= 1200 
~ 800 

° 400 

© 

> 0 

0 4 8 12 16 

Distance along beam path [cm] 



a) b) 

Fig. 3 Plots of a) CFD pressure data along laser LOS in vertical translation configuration (units are pascals) and b) temperature, pressure, and velocity 
data along LOS at the vertical center of duct. 


temperature, velocity, and pressure were extracted along the diagonal 
laser path indicated by the red arrows shown in Fig. 2b . 

The thin black box in Fig. 3a outlines the 1 mm region of the slice 
that was used to determine the CFD predicted conditions along the 
laser LOS in Fig. 3b. The variation of temperature, pressure, and axial 
velocity is plotted as a function of distance along the downstream- 
pointing beam for the sensor in the vertical translation configuration. 
These data were used with Eqs. (9-18) to compute the path- 
integrated lineshape for each LOS. The process was repeated at 
various vertical positions in order to compare with measurements 
taken using the TDLAS sensor. The path-integrated WMS lineshapes 
of Fig. 4 were calculated from the CFD data of Fig. 3b. The 1341 nm 
transition (£" = 1962.5 1 cnr 1 ) with a modulation index m = 0.9 is 
considered in this analysis. In the vertical translation configuration, 
the path-integrated absorbance is roughly 0.16. A significant fre- 
quency shift (Doppler shift) can be seen between the lineshapes 
produced on the downstream- and upstream-pointing beams (illus- 
trated in Figs. 2b and 2c). Note that the entire lineshape is positive, 
since the signal is the root sum square of the lock-in amplifier output 
[see Eqs. ( 12 ) and (15)]. 

The frequency shift can be measured at various points along the 
lineshape; however, to develop a high-speed, robust velocity detec- 
tion algorithm, only the data in the laser frequency scan from the 
central peak to the two adjacent valleys is considered. Ideally, the 
velocities measured on either the high- or low-frequency face of 
the lineshape would be identical. However, the presence of 
nonuniformities in the beam path can distort the lineshape, causing 
slight differences in these measured frequency shifts. 

Figure 5a shows the variation of velocity inferred from path- 
integrated lineshapes as a function of boundary-layer (BL) thickness, 
where the thickness of the boundary layer was defined as the distance 


to the point at which the velocity reached 99% of the freestream 
value. In this analysis the same pressure, temperature, and mole 
fraction gradients were used for all boundary-layer thicknesses. We 
will first consider the top two curves in Fig. 5a (BL = 0.9 cm thick), 
which are produced with the data in Fig. 3b. It is clear that the 
simulated path-integrated velocity measurement is lower than the 
core velocity shown in the lower plot of Fig. 3b (1600 m/s). This is to 
be expected, since the measurement is obtained from a path- 
integrated lineshape, produced as the laser beam travels through two 
low-velocity boundary layers and the core flow. Indeed, when the 
previously described measurement technique was applied to 
lineshapes simulated with no boundary layers, the core velocity was 
recovered at all measurement points along the lineshape. Intuitively, 



Fig. 4 Frequency-shifted path-integrated WMS-2//1/ lineshapes 
simulated from CFD data. Frequency shift corresponds to a 1600 m/s 
core flow. 
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Fig. 5 Plots of a) measured velocities from path-integrated lineshapes with varying boundary-layer thickness (the position along lineshape refers to the 
location on the lineshape used for Doppler shift measurement, with 0 corresponding to the valley and 1 corresponding to the central peak) and b) mean 
difference between measured velocity of Fig. 5a and core velocity versus combined boundary-layer thickness as percentage of an 18.7 cm path length. 


this suggests that the velocity inferred from LOS measurements 
should decrease as boundary-layer thickness increases, as illustrated 
in Fig. 5b. 

It should be noted that in addition to an overall decrease in 
measured velocity as the boundary-layer thickness increases, the 
difference in velocity measured at various points along the lineshape 
also increases (moving along abscissa of Fig. 5a). This is due to 
increased distortion in the lineshape as the low-velocity boundary 
layers become a greater percentage of the path length. This distortion 
also contributes to the increasing discrepancy between velocities 
measured on the high- and low-frequency faces as boundary-layer 
thickness increases. 

Figure 5b shows the mean deviation of the path-integrated 
measurement from the core velocity as a function of the combined 
boundary-layer thickness. The abscissa of Fig. 5b is computed as the 
percentage of the beam path along the 1 8.7 cm laser LOS that passes 
through the combined boundary layers. In the extreme case where 
roughly 40% of the laser LOS is contained in the boundary layers, 
only a 6% decrease in the path-integrated velocity measurement is 
expected. For the DCSCTF isolator, the CFD solution predicts that 
10% of the path length is contained in the boundary layers, and a 
decrease of less than 2% in the path-integrated velocity measurement 
is expected. It is clear that for the current transition and modulation 
index, the sensitivity of the TDLAS measurement to nonuniformities 
is quite low. Optimal line selection and optimal choice of modulation 
index can minimize the TDLAS sensitivity to nonuniformities; 
however, this discussion lies beyond the scope of the current paper. 

It is also seen that the decrease in the path-integrated velocity 
measurement is roughly linear with the boundary-layer thickness. 
This indicates that the path-integrated velocity measurement can be 
corrected in order to recover the core velocity. For example, with 
10% of the beam path passing through boundary layers, the simu- 
lated path-integrated velocity measurement is roughly 1.5% lower 
than the core value, and the TDLAS measurement can be adjusted 
accordingly. In the next section, this analysis of spatially resolved 
lineshapes simulated with the CFD solution is used to correct 
velocity measurements in the NASA Langley DCSCTF to within 
0.25% of the facility predicted value for the direct-coupled flow. 

III. Measurements in NASA Langley DCSCTF 

The DCSCTF at NASA Langley produces high-enthalpy flows for 
the testing of scramjet combustor concepts. H 2 -air combustion is 
used to heat the flow before expansion through a supersonic nozzle 


[24,1_5]. During testing of the TDLAS mass-flux sensor, a Mach 2.65 
nozzle was used to simulate Mach 6 and 7 flight enthalpies, 
producing the nominal test conditions summarized in Table L These 
conditions were computed using a 1-D thermodynamic equilibrium 
solver based on input mass fluxes to the combustor, nozzle areas, and 
heater pressures [16,f7]. The solver outputs predicted conditions at 
the nozzle exit based on the facility operating conditions for each run 
(input mass fluxes, heat losses, etc.). The DCSCTF flow tunnel is 
shown in a photograph in Fig. 6 with the major components labeled. 

The mass-flux sensor was mounted on a custom isolator section 
located between the nozzle and the fuel injectors in Fig. 6. Optical 
access via slot windows on the isolator section is shown in Fig. 7a. 
These windows cover 3 -mm- wide slots cut in the wall of the section 
to provide a crossing angle 26 = 90°, as illustrated by the arrows in 
Fig. 7a. 

All windows were 0.25-in. -thick BK-7 optical glass; window 
thickness was chosen such that internal reflections of the laser beam 
within the window would not be transmitted through the slots. The 
volume between the window and the slots was purged with heated 
nitrogen to minimize water vapor in the stagnant gases in the optical 
path through the isolator wall. The sensor optical components were 
attached to translation stages, as shown in Fig. 7b. The short optical 
paths external to the isolator were also purged with unheated 
nitrogen. The optical setup was translated on the vertical walls, then 
on the horizontal walls to enable measurements in both vertical and 
horizontal planes. 

A. Velocity Measured in the DCSCTF 

Figure 8a shows the time-resolved (50 Hz) velocity measured in 
the center of the channel with the sensor configured to probe the 
horizontal plane and a modulation index m = 0.9. After the initial 
startup transient of the tunnel (caused by a surge of flow after the 
hydrogen and oxygen supply valves fully open), the sensor measure- 
ment was in excellent agreement with the velocity computed from the 
facility predictive code mentioned earlier. The first plateau in velocity 
from 2.5 to 5 s corresponds to an H 2 and 0 2 mass-flow rate of roughly 
80% of the operating condition. The CFD solution for the isolator at 
the Mach 7 condition predicted the boundary-layer thickness to be 
roughly 7 mm. The uncorrected velocity measurements shown in 
Fig. 8a are roughly 1.7% lower than the facility predicted value; 
Fig. 8b shows the data of Fig. 8a with the boundary-layer correction 
and the improvement in agreement between measured and facility 
predicted velocities. 


Table 1 Conditions for Mach 6 and 7 set points at NASA Langley DCSCTF 


Flight enthalpy 

Static pressure, kPa 

Static temperature, K 

Water mole fraction 

Nozzle exit velocity, m/s 

Mach 6 

72 

714 

0.19 

1440 

Mach 7 

72 

935 

0.25 

1630 
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Fig. 6 Photograph of the DCSCTF at NASA Langley showing flowpath sections as labeled. 



Fig. 7 Mass-flux sensor installed on custom isolator section: a) sensor configured to probe vertical planes of the flowpath; arrows illustrate the beam 
paths and b) sensor configured to probe horizontal planes of the flowpath. 


By correcting for the presence of boundary layers, the difference 
between the sensor-measured velocity and the predicted velocity is 
reduced by more than a factor of 6, bringing the agreement to within 
0.25% in a 1630 m/s flow. Because of the necessity for accurate 
measurement of mass flux, this improvement in accuracy by 
accounting for nonuniformities along the LOS is quite significant. 
Measurements were also made with a modulation index m = 2.2 to 
investigate the improvement in the sensor performance from the 
optimization of the modulation index as described in [8], The use of 
an optimal modulation index reduced the standard deviation of the 
velocity after the startup transient by a factor of 3 (a = 14 m/s for 
m = 0.9 versus a = 47 m/s for m = 2.2). These data show that the 
velocity stability of flow in the tunnel and the sensor precision are 
better than 1% of the total 1630 m/s velocity after the startup 
transient. 


The sensor was then translated during a test across the channel for 
spatially resolved measurements; a scan of vertical planes from the 
left to the right sidewall of the channel (looking downstream) is 
shown in Fig. 9a and a scan of horizontal planes from top to bottom is 
shown in Fig. 9b. The sensor was first positioned within 0.5 mm of 
the wall at its starting position (top wall for vertical translation, left 
sidewall for horizontal translation) and measurements were taken at 
ten points as the sensor was translated across the full height or width 
of the channel during a single tunnel run. At each point in Fig. 9, the 
sensor measurements were averaged for 1 s, and the error bars shown 
are the 1 — a standard deviations of these average values. The total 
time for translation and data acquisition was 13 s. The optical 
alignment was degraded near the walls of the channel, where 
vibrations can block the edge of the laser beam; this reduced the 
signal-to-noise ratio of the transmitted intensity and produced a 




a) b) 

Fig. 8 Time-resolved velocity in the middle of the channel, horizontal plane: a) no correction and b) with correction for nonuniformity along LOS. 
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Fig. 9 Spatially resolved velocity (no correction applied) plotted from a) left to right of channel (facing downstream) in vertical planes and b) top to 
bottom of channel in horizontal planes. Solid data points indicate measurements taken during facility startup transient. 


corresponding increase in the observed standard deviation shown in 
the error bars. The spatially resolved CFD velocity data were path- 
integrated along the sensor beam path, and the resulting values are 
plotted in Fig. 9 for comparison. 

Overall there is good agreement between the CFD path-integrated 
solution and the measured values. The first three TDLAS measure- 
ments (solid data points in Figs. 9a and 9b) tend to be higher due to 
the startup transient seen in Fig. 8. In the vertical scan direction, the 
boundary-layer thickness predicted by the CFD (~9 mm) is 
confirmed by the TDLAS measurements. For the horizontal scan in 
Fig. 9a, the boundary-layer thickness may be underpredicted by the 
CFD simulation, although the TDLAS measurements agree closely 
in the core flow. These data demonstrate that the sensor has the 
capability to make precise, accurate, spatially resolved velocity 
measurements in the high-speed, high-temperature flow of the 
DCSCTF at NASA Langley. 

B. Temperature Measured in the DCSCTF 

Temperature was determined simultaneously with velocity during 
all runs. Both laser wavelengths were multiplexed onto the upstream 
and downstream beams, and temperature was determined by 
comparing the measured 1 /-normalized WMS-2/ peak values with 
simulations for both beam paths. The 1349 nm laser was modulated 
at m = 2.2 and the 1341 nm laser at m = 0.9 for velocimetry. As 
described in [8-11], the WMS-2/ signal is maximized at m = 2.2 
and hence changes slowly in response to perturbations such as flow 
nonuniformity; while only the 1349 nm laser was modulated at 
m = 2.2, this mitigated the influence of nonuniformity on the 
temperature measurement. Time-resolved temperature data are 
plotted in Fig. H) with the Mach 7 flight condition on the left and the 
Mach 6 flight condition on the right. These data were measured at the 


center of the flow for a horizontal plane at Mach 7 and a vertical plane 
for Mach 6. The temperatures for the downstream- and upstream- 
pointing beams are in excellent agreement, as expected. The 1 — <r 
standard deviation of the average temperature after the startup 
transient is ~10 K for both conditions. The gas temperature 
predicted by the facility model is also in excellent agreement with the 
measured data after the startup transient. More measurements are 
required to determine if the slow variation in gas temperature is a real 
facility effect or a measurement artifact, although the gradual rise in 
temperature during the run could reflect the change in wall losses as 
the test-section temperature increases. In Fig. 10b , the fluctuations in 
temperature are a result of the facility being unable to steadily 
maintain the Mach 6 flight condition (the nozzle used is designed for 
the Mach 7 condition). 

C. Mass Flux Measured in the DCSCTF 

The velocity and temperature data were combined with facility 
pressure measurements to determine the mass flux as the product of 
velocity and density. For these data the A.! ~ 1349 nm laser was 
modulated atm = 2.2andL 2 ~ 1341 nm was modulated at m = 0.9. 
The time-resolved mass flux is shown in Fig. IT, measured in the 
center horizontal plane for the Mach 7 flight condition in the left 
panel (Fig. 1 la ), and in the center vertical plane for the Mach 6 flight 
condition in the right panel (Fig. lib ). The data for mass flux from the 
downstream- and upstream-pointing beams are nearly identical, as 
expected from the excellent agreement in the temperature data from 
Fig. 10. There is also close agreement between the measured values 
and those predicted from the facility model; however, deviations 
occur near the beginning of the test due to the startup transient and 
toward the end of the test as the input gas mass flows decreased. The 
decrease in the measured mass flux was primarily caused by a 




a) b) 

Fig. 10 Gas temperature for downstream- and upstream-pointing beams a) in the center horizontal plane for the Mach 7 flight condition and b) in the 
center vertical plane for the Mach 6 flight condition on the right. Facility model value also shown. 
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Fig. 11 Mass flux using temperatures taken with downstream- and upstream-pointing beams: a) in the center horizontal plane for the Mach 7 flight 
condition and b) in the center vertical plane for the Mach 6 flight condition. Facility model value also shown. 
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Fig. 12 Spatially resolved mass flux at Mach 7 condition plotted from a) left to right of channel (facing downstream) in vertical planes and b) top to 
bottom of channel in horizontal planes. 


decrease in the 0 2 mass-flow rate supplied to the facility. The 
TDLAS measurements confirm that the mass flux was nearly 
constant during a run and demonstrate the sensor’s ability to detect 
slight deviations in mass flux caused by changes in the input mass- 
flow rates. 

Temperature and velocity data collected at various points within 
the duct were processed to produce the spatially resolved mass-flux 
measurements in Fig. 12 . The sensor measurements are in good 
agreement with the CFD solution; again, the three outliers at the far 
left of Fig. 12a and top of Fig. 12b are a result of increased mass flow 
during the startup transient. As mentioned previously, optical align- 
ment tends to be poorer near the edges of the slot windows, resulting 
in the higher standard deviations seen in the sensor measurements 
near the walls. Overall, the sensor measurements show that there is 
good spatial uniformity within the test section and prove that the CFD 
solution correctly models the DCSCTF flow conditions. 


IV. Conclusions 

TDLAS of water vapor at 1 .4 /xm was demonstrated to provide 
accurate, spatially and temporally resolved mass-flux measurements 
in a high-speed, high-temperature flow in the DCSCTF at NASA 
Langley. An analysis of the effects of nonuniformity in temperature, 
pressure, mole fraction, and velocity on the WMS-2//1/ lineshape 
was presented for the first time. By incorporating a CFD solution for 
the flow within the DCSCTF isolator section, path-integrated WMS- 
2//1/ lineshapes were produced, accounting for nonuniform 
conditions along the laser LOS. A correction to recover the core 
velocity from the path-integrated velocity as a function of boundary- 
layer thickness was developed. The sensor was then deployed in the 
DCSCTF at NASA Langley in a custom isolator flow section with 


optical access. Measurements of mass flux, temperature, and velocity 
were found to be in close agreement with both the facility predictive 
code (a 1 -D thermodynamic equilibrium solver) and a CFD solution. 
The TDLAS velocity measurement was corrected using the CFD 
path-integrated simulations to within 0.25% of the value predicted by 
the facility code. Finally, spatially resolved mass-flux measurements 
were shown to be in good agreement with the CFD solution, 
demonstrating the sensor’s capability to make accurate measure- 
ments with spatial resolution on the order of a beam diameter 
(1-2 mm). These results demonstrate the potential of TDLAS 
sensing for accurate mass-flux measurements in ground-test facilities 
and the potential for improving line-of-sight TDL measurements in 
nonuniform flowfields. 
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